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Abstract 

We calculate the mass of 0 ^ triple-valence-gluon resonance, the trigluon 
glueball, with QCD sum rules. Its mass is found to be approximately in the 
region between 1.9 GeV and 2.7 GeV with some theoretical uncertainties. 
Moreover, it is likely that the new BES measurement of the pp enhancement 
near threshold in the J/pj decays exhibits the behavior of this trigluon state. 
Our analyzes favor the baryonium-gluonium mixing picture for the BES ob¬ 
servation. 
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So far, one of the mysteries remaining in the Standard Model is the undiscovered 
elusive gluonium, or glueball, which is theoretically expected based on the non-Abelian 
and conhnement natures of Quantum Chromodynamics(QCD). Although glueball has 
been searched for by experiments for many years, there has been no clear evidence yet. Due 
to gluon self-interaction, the valence gluons inside the gluonium, suppose it exist, should 
be heavily ’’dressed” through the vacuum fluctuation. This kind of non-perturbative 
interaction, together with the bound state conhnement effects, makes the theoretical 
investigation pretty difficult. A even worse situation is that there might be a strong 
glueball-meson mixing, which may obscure the predicted glueball detection in experiment. 

SVZ QCD sum rule P has been applied successfully to many hadron phenomena, 
such as hadron spectrum and hadron decays. In this method, different interpolating 
currents are constructed corresponding to different hadrons, then the sum rule is estab¬ 
lished by constructing a correlation function and matching its operator product expansion 
(OPE) to its hadronic saturation. In the correlation function, the short- and long-distance 
strong interactions are separated by means of the OPE, where the former is perturbative 
calculated, whereas the latter are treated as non-perturbative fundamental parameters. 
Although this approach is plagued by its large theoretical uncertainties, it gives a model 
independent treatment of hadrons. At least, its results are qualitatively reliable and mean¬ 
ingful. Especially, before the lattice simulation, which starts from the QCD hrst principle, 
can make a more precise prediction on the physics of hadronic states, e.g. glueball, these 
kinds of studies are necessary. 

The two-valence-gluon resonances, the two-gluon glueballs(or bigluonium), have been 
studied extensively in the literatures. With QCD sum rules, the mass of scalar 
glueball 131311113011100113 im 113 and 2++ tensor glueball were calculated 0113 ) 
and their masses were predicted to be 1.6(±0.3) GeV and 1.7(±0.5) GeV, respectively. 
According to lattice simulations [T71 ITHl ITTH |23 HI]) the 0 glueball should be the third 
lightest pure gluonic state. Relevant theoretical investigations on this O"’*' state have been 
performed, but the studies are limited only to the two-valence-gluon glueball (the leading 
Fock state in Fock space expansion) j3 113113 ■ should be noted that the glueball 
can also be constructed by tree-valence gluons, which is different from the next-to-leading 
Fock state of the bigluonium in Fock space expansion. About this point, one can easily 
hnd similar cases in the quark hadron sector. 

In this work we will calculate the three-gluon O”’*' glueball mass by means of the QGD 
sum rules. The only independent triple gluon Lorentz irreducible and gauge invariant 
interpolating current for O”’*' glueball is found to be 




( 1 ) 


where. 



( 2 ) 
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fabc jg antisymmetric SU(3) structure constant, is the gluon field strength tensor, 
and three Gs are introduced here for the convenience of calculation. 







Figure 1: The typical Feynman diagrams of three-gluon glueball in the scheme of QCD 
sum rules, (a) diagram for the perturbative term; (b) for the two gluon condensate terms; 
(c) and (d) for the three-gluon condensate terms; (e) for the four-gluon condensate terms. 


The two-point correlation function for the sum rules is then constructed as: 

ntf) = *|dVe‘’-{o|r{i(i)j(o)}|o) , (3) 

where | 0 ) represents the physical vacuum. The Feynman diagrams of the correlation 
function in our calculation are shown in Figure 1. After a lengthy calculation, some 
formula are given in the appendix, the final result of the OPE for the correlation function 
up to order of four-gluon condensate is obtained 

n((5^) = ciqQ^ In ^ -|- boQ'^{asG‘^) + (cq -|-Ci In ^-^)Q^{gsG^) +(1^ \Osa) +do \Osb) , (4) 

where fi is the renormalization scale, > 0, and the coefficients are 


ao = — T^,bo = cq = —ISQttq;^, Ci = SdTro;^, = 0, (Iq^ = —9{47r)^as ■ (5) 

5 dTT 

The condensates are defined as 

{asG^) = , ( 6 ) 

{9sG^) = {9sr'^G;^Gl^G ;^), ( 7 ) 

{Osa) = , ( 8 ) 

(Osb) = {{asr’’^G';^Gl^r) . ( 9 ) 
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With the exact form of Eq. (ED, the dispersion relation with subtractions reads 

n(Q 2 ) = I*"” + (n(o) - Q^n'(o) + iQ‘‘n"(o) - lQ'=n"'(o)), (lo) 

where n(0), n'(0) and n"(0) are constants relevant to the correlation function at origin. 

On the other hand, in the narrow width approximation, the imaginary part of the 
correlation function can be saturated as: 

-lmn(s) = fM^6{s - M^) + p{s)e{s - Sq) , (n) 

TT 

where p(s) is the spectral function of excited states and continuum states above the 
continuum threshold sq, M is the mass of the resonance (glueball) and / is the coupling 
dehned by: 

(0|j(0)|G) = /M^ (12) 


To take control of the contribution of high order condensates in OPE and the con¬ 
tribution of high excited resonance and continuum states in the hadron side, the Borel 
transformation is performed 


B = 1- M ' 

" “ (n-1)! \dQ^^ 


to the sum rule, and the moment Rk is obtained 


Rk{T,so) = -B^ 


:-«")* (nw") - 




r L V TT Jso s'^{s + Q‘^) 

1 roo poo 

— / s^e“^'^Imn(s)(is — / p{s)ds . 

TT «/ 0 sq 

Employing the quark-hadron duality approximation 


( 13 ) 


( 14 ) 



1 

s^e-"Xs)ds ~ - / s^e-"OmnP'^^'(s)ds , 

7^ JSQ 


( 15 ) 


we have 

1 /■®o , 

Rk{'r,So) = — s^e ^'^Imn(s)(is . 

TT JO 


( 16 ) 


Then, the glueball mass can be extracted from the ratio of the above moments, i.e.. 


M^{so,t) 


Rk 

Rk -1 


(17) 
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The first two moments in Eq. m are 


Ro{so,r) = - ^[1-p4{soT)] + ^{gsG^) , 

ro 

_ / ^ 120 ciq p , , 2 c]^ . ^‘i\ 

Ri{so,r) = - ^[1-P5{sot)] +—{gsG^), 

where 


Pk(,x) 


k 


e-Y. 

n=0 



(18) 

(19) 


( 20 ) 


From what obtained above, we can easily get the analytic function of glueball mass relying 
on the Borel transformation parameter r and threshold Sq- 

For numerical analysis, the input parameters in our calculation are taken to be 


(asG^) = 0.06 GeFS 
(S.G®) = {0.27 GeV^-){a,G^-) , 

{0,1,} = ^{a.GY, 

16 

A.JPS ~ 300MeV , 

— 4:71 

llK^Afjs) ■ 


( 21 ) 


The three gluon condensate in dilute gas instanton model is used here for lack of 
direct extraction from experimental data. It is found that the variation of three gluon 
condensate changes the predicted glueball mass small. Other parameters are commonly 
used in literatures. The numerical dependence of glueball mass on the r is presented in 
Figure 2. There is not a exact stable window in the hgure. However, the glueball mass 
decreases slowly with the increase of Borel parameter r, an approximate glueball mass 
could be determined. The excited states and continuum threshold parameter sq = 6.5 
GeV^ is chosen according to the method in Ref. ng. To determine this sq, the suitable 
region for r lies between 0.2 GeV~^ and 0.7 GeV~^. At every r in this region, the predicted 
glueball mass varies slowly with sq. Therefore, the glueball mass could be approximately 
regarded as falling in the scope of 1.9 GeV to 2.7 GeV. 

The mass of 0 ^ three gluon glueball obtained here is not as large as one expected 
(the mass of O^’*' three gluon glueball is predicted around 3.1 GeVl^^l, where the nu¬ 
merical analysis is a little different from ours). There are several possible reasons for it. 
One possibility is that we have under-estimated the mass of three gluon glueball for the 
theoretical uncertainties. This under-estimation may originate from several ways. 

In the sum rules constructed with gluons current, two gluon condensate is expected to 
give a large contribution, which may change the mass determination. However, the leading 
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order two gluon condensate does not enter into the sum rule after Borel transformation 
both in our calculation and in Ref. [221 • If one wishes to keep two gluon condensate in 
SVZ sum rules and to see how large they will contribute, it is necessary to calculate the 
next leading order contribution of two gluon condensate or to know the behavior of the 
correlation function at origin. 

As pointed out in Refs. iniiisi, in the pseudoscalar glueball (two gluons) and rj' chan¬ 
nels, the topological charge screening contributions to the sum rules are found to be large. 
This screening effect generated from the interaction between the topological charges is ex¬ 
pected to exist also in the pseudoscalar three-gluon channel though its exact contribution 
has not yet been studied. 

In addition, in the pseudoscalar flavor singlet channel, the mixing among glueballs 
(three-gluon and two-gluon glueball) and normal mesons {rj and p') may change our result. 
All these contributions are beyond our scope in this paper. Certainly, the variations of r, 
Sq and other input parameters will result in some uncertainties also. 

Another possibility is that we have touched the reality. Our result for the O”'*' glueball 
is consistent with the recent lattice result (2.56 GeV)j2I]. As we know, in the lattice 
calculation, the predicted O”’*' glueball includes two gluon glueball, three gluon glueball 
and glueball with more than three gluons. If this is the truth, it implies that one additional 
dynamical gluon attributes less than one GeV to the glueball mass, not as what normally 
think in the hybrid case. To clarify this point of view, more investigations in other models 
are required. 



Figure 2: The trigluon glueball mass square dependence on the Borel parameter r. The 
excited states and continuum threshold parameter sq is taken to be 6.5 GeV^ 
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In recently, the BES Collaboration has observed an enhancement at the pp threshold 
with the quantnm nnmber O"’*' j2Sl and a resonance X(1835)[21] with undetermined. 
It is still unclear whether these two enhancements are the same. These observations have 
stirred many explanations such as the pp baryonium state 122], or the gluonic state pnicii 
[2B] • hr our understanding, the pp enhancement is probably a mixture of baryonium state 
with the glueball(both two- and three-gluon states), which is also proposed in Ref.j2Sj. 
However, we think that if the baryonium-gluonium mixing picture is the reality for the 
BES observation j23j. there should be large component of three-gluon glueball inside. The 
arguments for this idea are: 

1) The observed pp enhancement is close to the mass scope of O"’*' glueball in our 
calculation. 

2) The new observation exhibits large 7 -|- pp(and also three mesons) branching ratio 
relative to the two meson production processes in comparison with the decays of other 
known states. This unique character may hint that the new state has relatively strong 
coupling to six-quark hnal states Q. In the picture of the three gluon glueball, this is 
easy to be understood, since the three valence gluons can split into six light quarks easily 
and form pp hnal states. 

3) Theoretically, the rate oi J/p) radiative decays to three gluons is in a similar order 
of magnitude as to two gluons, while we should enforce a proper physical cut to avoid the 
infrared singularities to draw such a conclusion inoi. 

4) The missing AA threshold enhancement implies the obvious havor SU(3) breaking 
and favors the idea that what BES observed is a mixed state of pp baryonium and glueball, 
and the latter with a mass closing to the pp threshold. The pure baryonium interpretation 
meets the difficulty of large SU(3) breaking, while pure gluonic state interpretation has 
an overshoot mass. The three-gluon glueball-baryonium mixing(oscillation) picture is a 
simple and natural one to this understanding. 

In conclusion, in this paper we calculate the O”’*' three-gluon glueball mass by means 
of QCD sum rules. We hnd that its mass lies in the region of 1.9 ~ 2.7 GeV, which is 
consistent with the lattice simulation result. Since our calculation is performed in the 
leading order, the instanton and topological charge screening contributions have not been 
included in, the mixing effect among glueballs and normal mesons has not considered, 
and the input parameters remain some uncertainties, the mass prediction can not be very 
precise. In our leading order analysis, only the non-perturbatively triple-gluon condensate 
gives contributions to the glueball mass, whereas the two- and four-gluon condensates 
vanish after the Borel transformation. In further effort on the aim of making an accurate 
prediction, one should include the higher order corrections from both perturbative and 
non-perturbative sectors, in which the two- and four-gluon condensates contribute. The 
instanton and the topological charge screening contributions should be taken into account. 
Furthermore, the mixing effect among glueballs (two-gluon and three gluon glueball) and 
normal mesons (7 and p') should be considered. 
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In fact, what we calculated in this work is partially motivated by the new observation 
at BES. From the naive arguments listed above, we think the new BES observation of 
an enhancement at pp threshold might be a mixed state of baryonium and three-gluon 
glueball. Nevertheless, if our conjecture is true, how much the gluonic, or baryonium, 
content engages in the observed enhancement still needs more investigations, in both 
experiment and theory. 
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Appendix A 






d^pi d^p 2 d^ps {2Tr)^6^{q - Pi - p 2 - Ps) 


9 2 2 

P 1 P 2 P 3 


(27r)^ (27r)^ (27r)^ 

X{(-D^ + 8D- 14) [pKpi ■ p2f + pI{p2 ■ Psf + pUpi ■ Pzf 

+ (-D^ + 9L>2 - 29L> + 32)pIpIpI + (16 - QD){pi ■ P2)iP2 ■ P3)iPi -Ps)} 


3a3 /I 87 , ' 

10^^ (^ ^ 40 47ri/2 


n4(g^) = 

'' I_I I_I 

.3 (D - 2 ) 


= 27jg; 
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D(D-l) 

f d^pi 1 


{GG) 


{2ti)^ Pl^P2'^ 
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[{D^ -QD + 8)(pi ■ P 2 Y + {D^ - 7D^ + 17D - 14)pi^p2'] 


P2=q-pi 


13 


n 6 (gg = ^ SniQ^ 


(A.l) 


n=l 
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^l(g2) = lgig7j-abcjijkjamn 




= I 


.27 gUgsGGG) 


J L 


d^Pi 1 


2 D{D -1){D-2)J {27T)Dpfpl 
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P2=q-pi 


4^ + 1) 
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-«4g;{9.GGG> [ d°pi (D -l)(D- if Pi ■ (9 - pi) 
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D{D-1){D-2)J (27r)^ 
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PiPl 


p 2 =q-pi 
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J L 




= 0 


^5(g2) = ^Q-g7jabcjijkjamn 


X /<;'‘xe‘«-'[i£„g,4.1™(0)GS„p(x)l|i€„„„^?(0)G5,^(x)]{0|G;„(x)GS9,(0)GS,„(0)|0> 
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... , 6 /^ D-‘2 f d^pi Pi-{q- Pi) 

-t5AgAg,GGG)^-J 

^ \ 


-54TOfQ4i-7£-lnj^ + | 


V“'V Pl)^ 

^ + I) (9.GGG) 


5'6(Q2) = 


X fce*^-[-e«/3.^A™(0)G^,^(x)][Go^^jO)GS^^(x)](0|A-(0)G^o/3,(0)G^^^ 


= 0 


Sj{Q'^) = SGglr’^'^P'^f-^^jjd^xd^yG'^-^ 
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J L 


-i54gl{GGG) r d^p 2 1 


D{D-1){D-2)J {2ti)d plpi 
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7 rabc fijk rdst f f ^4, lA iq-x 


[2{2D - 5)piP2 + ‘^{D - 4)(pi ■ P 2 ) ] 


Pl=q-P 2 


SsiQ^) = 369.7“r*/"‘ / dXd'j/e*''- 


x|G;^„(9:)G‘„^(0)][GS„,(x).l”(9)][Gtj(!,)GS^(0)]{0|G5,^(x)GS,,„(0)A‘(9)|0> 


J L 

D, 


-27igl{GGG) r d^p 2 1 


D{D - 1){D -2) J {27i)^pjp^2 ^ 
= -27TialQ^{gsGGG) 


A{D - A){p 2 ■ pi) P 2 + 4(2D - 5)p^P2 


\Pl=Q-P2 


Sg{Q^) = 36glr’^J^^’^f^^JJd^xd^yG‘^-^ 

><[G^o,uix)&oaym[Gl,{x)Gt^s{y)M^^^ 


27igl{GGG) r d^p 2 1 


D{D - 1){D -2)J {2'K)Dplpl L 
= -27nalQ^{gsGGG) 

Sio{Q^) = ISglr’^J^^’^f^^JJd^xdSe^x^ 


4(4 - D){p 2 ■ pi) P 2 + 4(5 - 2D)p^P2 


pi=q-p2 


A&oaymM][Gl,{x)Gl^,{0)][Gl^^{x)^^^^^ 


= 0 
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Sn(O^) = fjd'‘xd-‘ye‘'<'‘ 

x|GL«(0)GL(!/)]|GS„p(i)GS^(0)]|G5,^(i)GS,J0)|{0|G;^„(i)yl‘(j/)yl<(s)|0> 

I_I I_I I_I 

= 0 

Su(Q^) = j[d‘xd‘ye‘'<-^ 

x|G;„„(i)/lS(!/)|[GS„^(i)GJ^,(0)|lGS„.(i)GS,„(0)](0|GJ„j(0)Gj„j(j/)yi;(s)|0> 

I_I I_I I_I 

= 0 

Si3(Q^) = 9glf‘^‘’f’\f‘‘” jjd-'xd*ye‘‘’=‘ 

x|G;„„(i)Gt,(j/)]|GS„,(i)Gj^,(0)]|GS„(i)G5;„(0)|{0|GJ„^(0)/lJ(j/)3l5(s)|0> 

I_I I_I I_I 

= 0 

neW") = -54TO’O"0-7E-lnj^ + 0 (s.GGG) 
n8(Q") = 9igtr'^r’‘ [d‘x e-''-'|GS„„(i)G‘„s(0)|{0|G‘,,(x)GS„(x)G;^,(0)G;;„(0)|0) 

■’ I_I 

= -Q{Anfas{{aJ-^'^Gl,GlX) 

Field strength tensor: 

GlAx) = Gl^Ax)+gsr^^A%x)Al{x) 

where 

GlfxAx) = dfxAKx) - d^A'^^ix) 

In coordinate gauge: 

Al(x) ~ ^x-Gl^iO); AliO) ~ 0 (A.2) 

G;„(i) = Gt^M + jg.r'‘‘xVG';^(0)GUO) 
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Some contractions: 


G%{0) = G^M = 


(A.3) 


GlA^)Gl^^{y) 

I_I 

where 


■ d‘^p —i5° 
(27r)"^ 


^VOL 


p{p)' 


,-ip-{x-y) 


y{p) = PyPa9uy + PvPygyia “ P^iPyP ipa PuPaQ ^0 
I _ I G") Pi 

P pA _ ^jm 

GW^(y) = J- p-,3,0)e-'’’'^-'''' 
GS,..(x)GU(!') = / 7|A^Ap.(p)e-''-''> 

I_I F 

where 


^y,ypij{p) PpPpdi'fT “ 1 “ PuPaQpp PpPaQup PvPpQpa P {dpadup PppSua) 

Some gluon condensates: 

r‘(0|G;,(0)G^(0)|0) = - 9 p„Jpp){GG> 

5'^(0|G“^(0)Gp^(0)|0) - - n ( 9 pp 9 p<^ ~ 9 pf 9 pp)(GG) 

/“’“(0|G;,(0)G^(0)G“„(0)|0) = p(g_i^)(j,_2) t?ppppp(GGG) 

where 


rjiS 

Xpupaap 


9ppQvaQap 9pp9i'P9cra 9pcr9iPa9pP A 9p<y9vP9pOi 

9pa9up9crP S“ 9pa9ua9pp "h 9pP9up9cra 9pP9vcr9pa 


Four-gluon condensate: 

r^ 7 ‘^'^^( 0 |G“,( 0 )Gj,( 0 )G^^( 0 )G:;,( 0 )| 0 ) 

9pp9i'a9cr'y9ps 9pp9va9a&9P'y 9pp9i'p9a'y9aS 9pp9i'P9crS9a'y 

9pp9iy^9cra9pS 9pp9iy^9a-p9aS “1“ 9pp9iy59cra9Pj 9pp9iy59a-p9a'y 

9pcr9ipa9p-f9ps A 9pcr9iPa9pS9P'f S“ 9p<j9iPp9pj9aS 9p<j9iPp9pS9a'y 
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~^9fJ.crQu'fQpaQfBS QpaQu’yQpl39a6 9p(T9i-'S9pa9f3'r “ 1 “ 9p<79v&9p99a'l 

9pa9Vp9CF'y99^ 9pa9up9crS99pa9i>cr9P'y99^ 9pa9iycr9ps99j 

~^9p99up9cr'y9<y-& 9p09i'p9cr59a'y 9pfi9lycr9p^9aS 9pp9i^cr9pS9a'y 

~^9p'y9iyp9cra9pS 9p'y9iyp9crp9ciS 9P'y9i>cr9pa9P& 9p'y9i'(j9pp9oi5 

9p&9i'p9(Toi9P"i “ 1 “ 9pS9vp9cFp9ct') 9p59i'cr9pa9p^ 9pS9ua9pp9ce"i} 

^[i9pa9uP 9pp9va){.9p^9cr5 9pS9cr'f) {,9p'y9i'S 9 pS9i''Y^i,9 pa9<7P 9pp9a-a}] 

where 

{D + 

{D + 2)D{D - 1){D - 2){D - 3) 

+ iD- 

{D + 2)D{D - 1){D - 2){D - 3) 

In 4-dimension: 
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